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Abstract. Tiie initial temperature T;, energy density Ei, and formation time Ti of the initial 
state of the QGP formed in the heavy ion collisions at RHIC and LHC energies are determined 
using the data driven Color String Percolation Model (CSPM). Multiparticle production by 
interacting strings stretched between projectile and target form a spanning cluster at the 
percolation threshold. The relativistic kinetic theory relation for rj/s is evaluated as a function of 
T and the mean free path (Xmfp) using data and CSPM. ri/s(Ti, \mfp) describes the transition 
from a strongly interacting QGP at T /Tc ~ 1 to a weakly coupled QGP at T /Tc > 6. We find 
that the reciprocal of 77/s is equal to the trace anomaly A = e — ZP/T'^ which also describes the 
transition. We couple this initial state of the QGP to a ID Bjorken expansion to determine the 
sound velocity of the QGP for 0.85 < T /Tc < 3. The bulk thermodynamic quantities and 
the equation of state are in excellent agreement with LQCD results. 



1. Clustering of Color Sources 

Multiparticle production is currently described in terms of color strings stretched between the 
projectile and the target, which decay into new strings and subsequently hadronize to produce 
observed hadrons. Color strings may be viewed as small areas in the transverse plane filled with 
color field created by colliding partons. With growing energy and size of the colliding system, 
the number of strings grows, and they start to overlap, forming clusters, in the transverse plane 
very much similar to disks in two dimensional percolation theory. At a certain critical density 
a macroscopic cluster appears that marks the percolation phase transition. This is the Color 
String Percolation Model (CSPM) [H [2]. The interaction between strings occurs when they 
overlap and the general result, due to the SU(3) random summation of charges, is a reduction 
in multiplicity and an increase in the string tension hence increase in the average transverse 
momentum squared, {pD- We assume that a cluster of n strings that occupies an area of 5„ 
behaves as a single color source with a higher color field Qn corresponding to the vectorial sum 
of the color charges of each individual string Qi. The resulting color field covers the area of 
the cluster. As Qn = X^i Qi; f^e individual string colors may be oriented in an arbitrary 
manner respective to each other , the average QuQij is zero, and Q"^ = nQ\. 

Knowing the color charge Qn one can obtain the multiplicity jj, and the mean transverse 
momentum squared {pi) of the particles produced by a cluster of n strings [2] 




(1) 



where //q and {pf)i are the mean multiphcity and (p^) of particles produced from a single string 
with a transverse area Si = vrrQ. This implies a simple relation between the multiplicity and 
transverse momentum fJ'n{Pt)n = ''^/^o(Pt)i) which means conservation of the total transverse 
momentum produced. 

In the thermodynamic limit, one obtains an analytic expression [U [2] 

where F^) is the color suppression factor and ^ = is the percolation density parameter. 

Eq.(l) can be written as ;U„ = n-F(^)//o and (pj)n = {pi)i/P{0- The critical cluster which 
spans Sn, appears for ^ 1-2 [3]. It is worth noting that CSPM is a saturation model similar 
to the Color Glass Condensate (CGC), where {Pt)i/F{0 plays the same role as the saturation 
momentum scale in the CGC model [H[5]- 

2. Experimental Determination of the Color Suppression Factor F{^) 
The suppression factor is determined by comparing the pp and A+A transverse momentum 
spectra. To evaluate the value of from data for Au+Au collisions, a parameterization of pp 
events at 200 GeV is used to compute the pt distribution 

dNjdp^t=a/{po+Ptr (3) 

where a is the normalization factor. pQ and a are parameters used to fit the data. This 
parameterization also can be used for nucleus-nucleus collisions to take into account the 
interactions of the strings [2] 

dNjdp^, = , ^ (4) 

The color suppression factor F{S^) is related to the percolation density parameter 



i^(0 = y^— (5) 

3. Energy Density e 

The initial energy density above Tc is given by [6] 

^ '-'n'pro 

To evaluate we use the charged pion multiplicity dNc/dy at midrapidity and 5„ values 
from STAR for 0-10% central Au-Au collisions with ^/snn =200 GeV We can calculate 
(pt) using the CSPM thermal distribution Eqs.(9) and (10). For 0.2 < pt < 1.5, (pt) = 
0.394 lb O.OOSGeV, adding the extra energy required for the rest mass of pions at hadronization 
{mt) = 0.42 lb O.OOSGeV. The error on (pt) is due to the error on Tj [8]. 

The dynamics of massless particle production has been studied in QE2 quantum 
electrodynamics. QE2 can be scaled from electrodynamics to quantum chromodynamics using 
the ratio of the coupling constants [9]. The production time Tpro for a boson (gluon) is [10] 

_ 2.405;i 



For Au-Au collisions at ./IWn =200 GeV Eqs. (6) and (7) gives £i = 2.27±0.16 GeV/Zm^ at 
^=2.88. In CSPM the energy density e is proportional to S^. From the measured value of ^ and 
e it is found that e is proportional to for the range 1.2 < ^ < 10.56, = 0.788 ^ GeV/ fm^ 
[71 lllj as shown in Fig.l. 
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Figure 1. Energy density e as a function of the percolation density parameter ^. The value for 
LHC energy is shown as blue square [3 Hi] . 



4. Determination of the Temperature 

The connection between the measured ^ and the temperature T(^) involves the Schwinger 
mechanism (SM) for particle production. The Schwinger distribution for massless particles 
is expressed in terms of [lOl |9] 

dn/dpf ~ e-"P?/^' (8) 

where the average value of the string tension is (x^). The tension of the macroscopic cluster 
fluctuates around its mean value because the chromo-electric field is not constant. The origin of 
the string fluctuation is related to the stochastic picture of the QCD vacuum. Since the average 
value of the color field strength must vanish, it can not be constant but changes randomly from 
point to point [12]. Such fiuctuations lead to a Gaussian distribution of the string tension for 
the cluster, which transforms SM into the thermal distribution |jl2j 

dn/dpl ^e'"'^^^ (9) 
with (x^) = 7r(p^)i/F(^). The temperature is expressed as [5| [13] 

no = (10) 

Recently, it has been suggested that fast thermalization in heavy ion collisions can occur through 
the existence of an event horizon caused by a rapid de-acceleration of the colliding nuclei |14j. 
The thermalization in this case is due to the Hawking-Unruh effect |151ll6j. In CSPM the strong 
color field inside the large cluster produces de-acceleration of the primary qq pair which can be 
seen as a thermal temperature by means of the Hawking-Unruh effect. The string percolation 
density parameter ^ which characterizes the percolation clusters measures the initial temperature 
of the system. Since this cluster covers most of the interaction area, this temperature becomes 
a global temperature determined by the string density. In this way at = 1-2 the connectivity 
percolation transition at T(^c) models the thermal deconfinement transition. 



We adopt the point of view that the experimentally determined universal chemical freeze-out 
temperature (Tj) is a good measure of the phase transition temperature, Tc. {pf)i is evaluated 

using Eq.(lO) at = 1-2 with Tj = 167.7 ± 2.6 MeV ^7\. This gives ^J{p^)l = 207.2 ± 3.3 
MeV. This calibrates the CSPM temperature scale. 

Figure 2 shows a plot of e/T'^ as a function of T/Tc. The lattice QCD results are from 
HotQCD Collaboration [18]. It is observed that at LHC energy the CSPM results are in 
excellent agreement with the lattice QCD results. The lattice and CSPM results are available 
for T/Tc < 2. Table I gives the CSPM values ^, T, e and rj/s at T/T^ = 0.88, 1, 1.16 and 1.57. 
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Figure 2. The energy density from CSPM versus T/T^^^^^ (red circles) and Lattice QCD 
energy density vs T/T^^'-'^ (blue dash line) for 2+1 flavor and p4 action [T8] . 



5. The Shear Viscosity to Entropy Density ratio r]/s 

The relativistic kinetic theory relation for the shear viscosity over entropy density ratio, rj/s is 
given bv [T91I20] 

V ^ TXmfp 

s 5 

where T is the temperature and Xmfp is the mean free path given by 
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^"^^^ ^ (i) ^^^^ 

n is the number density of an ideal gas of quarks and gluons and atr the transport cross section 
for these constituents. 

After the cluster is formed it behaves like a free gas of constituents. Eq. (10) can be applied 
to obtain the shear viscosity. In CSPM the number density is given by the effective number of 
sources per unit volume j21j 

^sources /-, o\ 

L is the longitudinal extension of the source, L = 1 fm [13]. The area occupied by the strings is 
related to ^ through the relation (1 — e~^)Si^. Thus the effective no. of sources is given by the 



Table 1. The measured percolation density parameter temperature T, energy density e and 
■q/s for the meson gas [22j; the hadron to QGP transition; Au+Au at 200 GeV and Pb+Pb at 
2.76 TeV . Au+Au is for 0-10% and Pb+Pb is for 0-5% central events. 



System 




T (MeV) 


e{GeV/fm-') 


r?/s 


Meson Gas 


0.22 


150.0 




0.76 


Hadron to QGP 


1.2 


167.7 +2.6 


0.94+0.07 


0.240+0.012 


Au+Au 


2.88+0.09 


193.6 +3.0 


2.27+0.16 


0.204+0.020 


Pb+Pb 


10.56+1.05 


262.2+13.0 


8.32+0.83 


0.260+0.026 



total area occupied by the strings divided by the effective area of the string SiF{^). 

N - - (14) 

sources — SiF{^) \ ' 

In general Ngources is smaller than the number of single strings. The number density of sources 
from Eqs. (12) and (13) becomes 

In CSPM the transport cross section atr is the transverse area of the effective string SiF{^). 
Thus atr is directly proportional to F{^) and hence to t^. The mean free path is given by 

L 

= (1 _ (1^^ 

For a large value of ^ the Xmfp reaches a constant value, rj/s is obtained from ^ and the 
temperature 

^ = (17) 

s 5(1-6-5) ^ ' 

The behavior of ry/s is dominated by the fractional area covered by strings for £, < ^c- This is 
not surprising because r//s is the ability to transport momenta at large distances and that has 
to do with the density of voids in the matter. 

Figure 3 shows a plot of ry/s as a function of T/Tc. The estimated value of rj/s for Pb+Pb 
is shown in Fig. 3 at T/Tc = 1.57. The lower bound shown in Fig. 3 is given by AdS/CFT 
|23j . These results from STAR and ALICE data show that the -q/s value is 2.5 and 3.3 times 
the KSS bound [23]. 

The theoretical estimates of 77/5 has been obtained as a function of T/Tc for both the weakly 
(wQGP) and strongly (sQGP) coupled QCD plasma are shown in Fig. 3 |19]. It is seen that at 
the RHIC top energy r]/ s is close to the sQGP. Even at the LHC energy it follows the trend of 
the sQGP. By extrapolating the r//s CSPM values to higher temperatures it is clear that 'q/s 
could approach the weak coupling limit near T/Tc ~ 6. The CSPM rj/s value for the hadron 
gas is in agreement with the calculated value using measured elastic cross sections for a gas of 
pions and kaons |22j . 



6. The trace Anomaly A and the reciprocal of -q/s 

The A and rj/s both describe the transition from a strongly coupled QGP to a weakly coupled 
QGP. We find that the reciprocal of rj/ s is in quantitative agreement with e — SP/T"^ the trace 
anomaly over this wide temperature range. This result is shown in Fig. 4. The minimum in rj/s 
=0.20 at T/Tc 1-15 determines the peak of the interaction measure ~ 5 in agreement with the 
recent HotQCD values p^. 
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Figure 3. rj/s as a function of T/Tc. Au+Au at 200 GeV for 0-10% centrality is shown as solid 
black square. wQGP and sQGP values are shown as dotted blue and green lines respectively 
[19j. The red dotted line represents the extrapolation to higher temperatures from the CSPM. 
The hadron gas value for tj/s ~ 0.7 is shown as solid black circle at T/Tc ~0.88 [22^. 
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Figure 4. The trace anomaly A = e — 3P/T^ vs temperature. 



7. The Equation of State 

We use CSPM coupled to a ID Bjorken expansion. The input parameters are the initial 
temperature T, the initial energy density e, and the trace anomaly A are determined by data. 
The Bjorken ID expansion gives the sound velocity 

1^ = _C2/r (18) 
dT dT de 

— = (19) 

dr de dr 

de „ , , , 

- = -Tsir (20) 



where e is the energy density, s the entropy density, r the proper time, and the sound velocity. 
Since s = £ + P/T and P = {e - AT^)/3 one gets 

f»^Cj (21) 

From above equations can be expressed in terms of ^ 

Cl = (-0.33) - l) + 0.0191(A/3) [tT^, " JZ^) (^2) 
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Figure 5. The speed of sound from CSPM versus T/T^^^^^ {red circles) and Lattice QCD-p4 
speed of sound versus T/T^^'3'-^^(blue dash line)[ 



8. Conclusions and open questions 

Two central objectives in the experimental characterization of the QGP are the Equation Of 
State (EOS) and the shear viscosity to entropy ratio r]/s. We have found that s/tj is equal to 
the trace anomaly A = e — 3P/T^. We determine the bulk thermodynamics, e/T^, s/T'^ and 
the EOS (?g of the QGP as a function of r/rc(CSPM). The results are in excellent agreement 
with LQCD numerical simulations as a function of T/Tc(LQCD). CSPM predicts that the QGP 
will be formed in high multiplicity p-p collisions at LHC. The CSPM r]/s ratio can be used to 
test the quasi-particle hypothesis where the jet transport parameter is g = f T^^, which can be 
compared with q measured in jet quenching using the the ALICE EMCAL detector at LHC. 
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